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ABSTRACT 

We have shown that there exists low-frequency growing modes driven by a global temperature 
gradient in electron and ion plasmas, by linear perturbation analysis within the frame work 
of plasma Kinetic theory. The driving force of the instability is the local deviation of the 
distribution function from the Maxwell-Boltzmann due to global temperature gradient. 
Application to the intracluster medium shows that scattering of the particles due to waves 
excited by the instability is possible to reduce mean free paths of electron and ion down to five 
to seven order of magnitude than the mean free paths due to Coulomb collisions. This may 
provide a hint to explain why hot and cool gas can co-exist in the intracluster medium in spite 
of the very short evaporation time scale due to thermal conduction if the conductivity is the 
classical Spitzer value. Our results suggest that the realization of the global thermal equilibrium 
is postponed by the local instability which is induced for quicker realization of local thermal 
equilibrium state in plasmas. The instability provides a new possibility to create and grow 
cosmic magnetic fields without any seed magnetic field. 

Subject headings: cooling flows — instabilities — plasmas — galaxies: clusters: general 

1. Introduction 

One of the most important results obtained by Japanese/US X-ray satellite ASCA (Tanaka, Inoue, 
& Holt 1994) is the first spectro-scopical confirmation of existence of at least two different temperature 
X-ray emitting hot gas in the cluster central region (Fukazawa et al. 1994; Ikebe 1995; Ikebe et al. 1999); 
the cooler is about IkeV and the hotter is about 2keV-10keV. Although the interpretation of the result 
is still controversial (Fabian 1994; Ikebe et al. 1999), there is a common feature among the different 
interpretations, that is, a cool region with temperature Ti is surrounded by a hot region with temperature 
T2 maintaining a pressure equilibrium (see §6.2). Since the total thermal energy in the hot region is much 
higher than that in the cool region, the hot region can be treated as a heat reservoir. Suppose that the 
interface separating two plasma regions with different temperatures has a thickness of L and temperature 
variation in the interface region is almost linear in spatial scale. Assume that the electron mean free path Ae 
is smaller than L. The time scale t^v of the evaporation due to heat conduction can be roughly estimated 
by a time scale for electrons in the hotter region to diffuse across the interface, and is given as 



Tev ~ (Ae/wth,e) X (L/A,)^ = e^\L/vtk,e) ^ 10^ (i/4kpc) (A:Bre/4keV) ^^'^ (e/ 1 .0) " V^, 



(1) 
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where Uth,e = V^^BTeTwT is the electron thermal velocity, is the electron temperature, and rUe is the 
electron mass; e is defined by e = Ag/L. In the cluster central region, the electron mean free path due to 
Coulomb collision is obtained as Ac ~ 0.5(A:BTe/4keV)^(ne/10^^cm^'^)^^kpc (Sarazin 1988). Therefore, 
the time scale of the evaporation is about 10^ yr and is much shorter than the age of cluster of galaxies, 
which is about the Hubble time. For maintaining the interface at least as long as the age of the cluster, 
the reduction of the electron mean free path in four order of magnitude from the Coulomb mean free path 
to realize e of ^ 10~^ is required. Soker & Sarazin (1990) and Makishima (1997) proposed that the cool 
and hot phases are thermally insulated from each other by magnetic fields. However, their models involve 
another problem how such special magnetic field structures are obtained in the cluster central region. In 
this paper, we show that there is an alternative possibility to reduce electron mean free path in many order 
of magnitude without any special assumption on magnetic field structure. 

Ramani & Laval (1978) found a new plasma instability which may relate to the reduction of the electron 
mean free path, henceforce the reduction of the heat conduction. They showed that the temperature 

gradient leads to an anisotropic electron velocity distribution function and that the anisotropy of the 
velocity distribution function drives the instability like the Weibel instability (Weibel 1959; Fried 1959; 
Melrose 1986). They proposed that the chaotic magnetic field and electric field, produced by the plasma 
wave induced by the instability, scatter electrons, so that the electron mean free path is reduced in many 
order of magnitude. However, the instability found by them cannot be applied to reduce the electron mean 
free path in many order of magnitude because of the following reasons. Since they assumed that only 
electrons respond to the mode and ions were treated as fixed back ground particles, the phase velocity of 
the wave must be faster than the ion sound velocity Vth,i — y^^koTiJrni, where T^, n^, and m, are the ion 
temperature, ion number density, and ion mass, respectively. In the astrophysical plasma, the ion mass can 
be safely replaced by the proton mass. The unstable mode found by Ramani & Laval (1978) has non zero 
real part of the wave frequency with a phase velocity of evth,e- Therefore, the application limit of their 
analysis sets a relatively high lower limit on e as e > ^ra^jmi ^ 0.025 when Ti = T,. (hereafter we refer to 
this lower limit as the wall of the square root of the mass ratio). To explain the two-phase nature of the 
hot gas in cluster central region, the instability found by them is practically not useful and the mechanism 
which can break the wall of the square root of the mass ratio is required. 

The first application of the Ramani & Laval (1978) type instability to the astrophysical plasma was 

made by Levinson & Eichler (1992). They extend the Ramani & Laval (1978)'s analysis to include the 
non-zero background magnetic field and found the low-frequency unstable mode, similar to the Ramani 
& Laval (1978) type modes. Pistnner, Levinson & Eichler (199G) first applied the Levinson & Eichler 
(1992) 's results to the cluster cooling flow. However, they applied their results to reduce electron mean free 
path down to 10~^ and far less than the lower limit, bounded by the application limit of the wall of the 
square root of the mass ratio. Further they neglected the cluster gravitational force in their Boltzmann 
equation in spite of taking into account a pressure gradient with a scale height comparable to the cluster 
gravitational scale height. As a result, they implicitly assumed the existence of the steady electric field for 
the electric force acting on an electron to balance with the pressure gradient force acting on the electron. 
Since the pressure gradient force acting on an ion balances with the cluster gravitational force acting on 
the ion, the electric force acting on an ion is as strong as the cluster gravitational force acting on the ion. 
The existence of such a strong electric field in cluster of galaxies is unlikely. Further, the model adopted 
as the temperature distribution in the cluster central region obtained by White & Sarazin (1987) only 
describes the temperature distribution of the hottest phase gas. Therefore, the model adopted by them is 
inappropriate to answer why the multi- phase gas clumps can co-exist in the cluster of galaxies. 
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In this paper, we extend the Ramani & Laval (1978) 's analysis including the response of the ion to 
examine whether we can break the wall of the square root of the mass ratio. We limit our attention to the 
pressure equilibrium plasma; in other words, the scale of the interface is much smaller than the pressure 
scale height. Then, the gravitational force can be neglected. The plan of the paper is as follows. The 
readers who arc not interested in the details of the derivation of the dispersion relation can skip §§3 and 
4. In §2, we set the model, and summarize the basic assumptions and the application limit coming from 
the assumptions. In §3, we deduce the solution of the Boltzmann equation under the situations set in §2 
using the Chapman-Enskog expansion of the electron and ion velocity distribution functions. We deduce 
the dispersion relations imdcr low-frequency conditions in §4. In §5, we summarizes the characteristics of 
the instability and modes found in §4. The application of the results to the hot plasma in the central region 
of cluster of galaxies is discussed in §6. Finally, §7 summarizes another important implications from our 
results and remained problems. 



2. The model and application limit 

Consider the situation described in Fig. 1. Two plasma regions with different temperatures contact 
maintaining pressure equilibrium. The electron temperature in the transition region with a depth of 
L varies almost linearly from Ti to T2 (Ti < T2) along the a;-direction: VT{x)/T{x) = 6tL~^ where 
5t = {T2 — Ti)/T. In clusters of galaxies, the pressure scale height is about the core radius of the cluster 
mass distribution, say ~ 50kpc (e.g., Hattori, Kneib, & Makino 1999). Therefore, if the size of the interface 
is smaller than ~ lOkpc, the pressure equilibrium condition for the interface region can be safely applied. 
The ion temperature is assumed to be the same as electron temperature,Ti = Tg. The mean free paths Ae 
and Xi of electron and ion, respectively, are assumed to be defined by the scattering due to plasma waves. If 
the scattering of the particle due to electric field E of the wave is the dominant source, the mean free path 
of each particle could be described as follows. Assume that the electric field associated with plasma wave 
is coherent in the scale of order of the wave length A and is chaotic over the scale. Then the momentum 
change of the particles can be described by the random walk manner in the momentum space. The mean 
free path is defined as the length scale that the particles gain or loss its momentum comparable to their 
original momentum during they move the scale. The momentum change of electron during one step is 
TOgAue ^ ei?A/uth,e- The number of steps required for electron to gain or loss their momentum comparable 
to its original momentum is [meVth,e/{eEX/vth,e)]^ ~ {2kBTe/eEX)'^ . Then we obtain 



V eEX J \ eE 

Similarly, the mean free path for ion can be obtained as 



Since we are now assuming that electron and ion have the same temperature, the electron and ion mean 
free paths due to the scattering by the wave electric field are identical. When the chaotic magnetic field 
associated with the plasma wave is the dominant source of the particle scattering, the particle mean free 
path would be defined in the following way. The magnetic moment of the particles, /i = {l/2mev\)/ is a 
conserved variable in a magnetized plasma, where v±_ is the particle velocity perpendicular to the magnetic 
field, written as mev\ ~ 2kBTe, just after the perturbed magnetic field is induced. As the magnetic- field 
strength gets stronger, v±_ increases to conserve /U. The maximum allowable value of v±_ is ^J^k^f^Jm^^ 
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since the kinetic energy of the particle must be conserved. Then it follows that an electron at the place with 
a magnetic field strength of Bi is scattered back at the place with a magnetic field strength of Bcr ~ 3/2i?i 
owing to the mirror effect. The mean free path is then defined as the length scale along the magnetic field 
from the place where the magnetic field strength is Bi to the place where the magnetic field strength gets 
Bcr- Since the critical value B^ of the magnetic field strength, in other words the place of the bouncing, 
docs not depend on the particle species and only depends on their initial position and the magnetic-field 
configuration. Therefore, the mean free paths of ion and electron must be also the same in this case. These 
arguments show that the mean free paths of electron and ion due to scattering by plasma waves may have 
the same order of magnitude. Therefore, we assume that electrons and ions have a similar value of e. 

In this paper, perturbations with a short-wavelength wave, such that 

A < Ae, (4) 

are treated. Therefore, the collision term in the Boltzmann equation for the perturbed variables can be 
neglected. 

Self-consistency check of the obtained results with all these conditions constrain the allowable range of 
e. In addition, two further conditions provide important limits on e. First of all, the electron mean free 
path due to Coulomb collision must be shorter than L at least before the onset of the instability, otherwise 
there is no chance for the electron distribution function to have anisotropy as described in the next section, 
and hence no chance to have instability. This condition yields an lower limit on L as 

L > 0.5(fcB7;/4keV)2(ne/10-^cm-3)-ikpc. (5) 

Secondly, the growth time scale of the unstable mode, must be shorter than the time scale of thermal 
difl[usion when waves do not exist, in order for the waves to grow before the temperature difference will be 
erased owing to quick thermal conduction. This condition provides 

W;"^ < e~'^L/vth,e, (6) 

where only the case that the electron mean free path is shorter than L is considered. A further condition 
has to be required when applied to the clusters, that is the growth time scale must be shorter than the age 
of the object, conservatively the age of the universe: 

ojr' < 1/Ho = lO^O/i-Vr, (7) 
where Hq = 100/ikm/sec/Mpc is the Hubble constant. 



3. Distribution Function 

When e < 1, the discussion given in the previous section ensures that the solution of the Boltzmann 
equation of electron and ion can be expanded in powers of e (Chapman & Cowling 1960) 

fe = /xn,e + /i^) + /i') + ..., (8) 
fi = /m,i + /f^+/r+..., (9) 

where /m.e and /m.i are Maxwellian distribution functions for electron and ion, respectively; /i*^^ 
and f-''^ {k = 1,2, •••) describe the deviation of distribution functions from the Maxwellian in 
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order of e*^ for electron and ion, respectively. Since we assume that Ti{x) = Te{x) = T{x) and the 
charge neutrality must be maintained, the Maxwellian parts for both electron and ion are written as 
/m,e = no(a;)[7r?;th.,.(.T)]-3/2exp{-[?;/?;th,e(a;)]^} and /m,, = rip {x) [nvti^^i {x) ] '^/"^ exp{- [^ /t'th.i jx)]'^}, 
respectively. Il(>r(> is the electron number density, Vih^e{x) = ^J2k^T{x)Jrnl, and 

vth,i{x) = ^J2k■B.T{x)/mi. 



The pressure equilibrium assumption gives 

Vno VT 1 



5t- (10) 



np T L 

Once this condition is satisfied, the pressure is time independent everywhere, even if a secular variation of 
the temperature due to thermal conduction is taken into account. Then the time dependence of the density 
can be related to that of the temperature as 

no dt T dt' ^ ' 

The Boltzmann equation under the pressure equilibrium condition without background electric and 
magnetic fields leads to 

dfe^dfe fdfA 



dt dx \ dt ^ 

where {df/dt)c is the collision term and the equilibrium electric field does not appear because of the 
pressiirc equilibrium condition. The secular time variation of the distribution function is ascribed to the 
secular variation of the temperature due to thermal conduction. Since the cool region with a temperature of 
Ti is considered to be immersed in a heat bath with a temperature of T2, the temperature of the interface 
increases monotonically owing to the heat conduction. Hence, the time evolution of temperature can be 
described as 



Then, 



8t ^''^ 



fe - eSr^U (14) 



dfe 


1 


dT 


dt 




dt 



The simplest choice of the collision term is the Krook operator, such that 

'df, 
dt 



= -Mfe - /m,e), (15) 



where Ve = Vth,e/^e is the collision frequency assumed to be constant. Substituting the collision term by 
the Krook operator, ordering the Boltzmann equation in and provides 



Vp dx 



(16) 



It shows that the system can be treated as steady state in the first order of e but the secular variation 
cannot be negligible in the second order of e. The time derivative of the Maxwell-Boltzmann distribution 
function can be described as 
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The first- and the second-order solutions for the electron distribution function arc thus obtained as 
follows: 



Vth,e 



nieV 



fn 



f(2) _ 



2v' 



+ 



2v' 



(19) 
(20) 



where we have neglected the second derivative of T(x) with respect to x in derivation of the second-order 
solution. 

These equations show that the anisotropy in velocity distribution functions are induced by the existence 
of the temperature gradient. The first-order solution is an odd function of the velocity. On the other hand, 
the second-order solution is even in velocity. It is straight forward to make sure that the zero background 
electric field assumption is consistent with the zero electric current condition under the pressure equilibrium 
plasma condition, by checking (fev) — 0, where ( ) denotes the average over the velocity. 

By substituting nie to mi and Ve to Ui = Vi/Xt in equations ( p^ ) and (pO|), the first and the second 
order distribution functions fj;^^ and /^'^^ for ion are obtained. 



4. Dispersion relation 

In the short- wavelength limit, the coUisionless Boltzmann equation is a good representative for the 
evolution of the perturbed distribution functions, and the plain wave exp[i(k ■ f— Lut)] is a good representative 
for the perturbed quantities (e.g., electron distribution function as 6fe = fj: ^ exp[i(fc • r — wt)], electric 
field as E ~ exp[i(fc • r — ^^t)], and magnetic field as i? = exp[j(/c • r ~ t^t)]). From the linearized 
Boltzmann equation in terms of the perturbed quantities of both electron and ion, the equations for the 
perturbed distribution functions of electron and ion are obtained as 

-t{uj-k- v)fj:^^ 
-i{uj - k ■ v)f^^^ 

As shown below, the low-frequency condition requires that the charge density be small, implying 
that the electric current has to be almost transverse for purely electronic oscillations. Therefore, only 
two modes are expected. Without loss of generality, the wave vector k can be taken in the x-y plane as 
kx = kcos9, ky = fcsin6', and fcz — 0. We expect one mode with the magnetic field along the z-direction 
(mode (1): ^ By ~ and Bz 7^ 0), and the other one with the magnetic field in the x-y plane (mode (2): 
Bx, By 7^ and B^ — 0). We thus denote the component of a vector along the direction of the wave vector 
by a fc-subscript and the component along the perpendicular direction of the wave vector in the x-y plain 
by a _L-subscript. The x-component of velocity, for example, can be expressed as Vx = Vk cosO — v± sinO in 
terms of Vk and v± . 

The Faraday's law in Maxwell's equations leads to -E^ = (i.e., E^ = {Ek, E±, 0)) and kE± = {uj/c)Bz 
for mode (1). For mode (2), div_B = yields Bk = (i.e., B^ — (0, i3j_,0)), and the Faraday's law and 



— {Ej: + vx ^ 
nif, \ c 



— [Er + vx ^ 



dfe 
dv 

dh 
dv 



(21) 
(22) 
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B^ = lead to E± = and kE^ = 
for mode (1), and 

{k-v~uj)ff:,^ 

for mode (2). 



-(lu/c)Bj_. Then we obtain 



+ - 



le 
me 
ie 
TO,: 



c 



dvk 



El - 



c 



ie 

TO,; 



c 



5ie 

dv I 



5^ 



(23) 
(24) 

(25) 
(26) 



To first order in [j/(fctith,e) and w/(fcwth,i), (fc-w-w)"^ = {l/k)P{l/vk) + {uj /k'^)P{l/vl) + i{T: /k)5{vk)- 
i{'K/k){uj/k)(dS{vk)/dvk), where P denotes the principal value and the signs in front of the delta functions 
reflect the causality condition. Then, with the help of the Faraday's law, equations (|2^) and (p3) yield 



le 

TOe 



(Ek + —B, 
V c 



2vk , eSj 



2'yth,, 







cos 6 









2vk 



sin 61 5 



2v^ 



"th,e , 



„ ^ (wfe cos 6 — v± sin ( 



2 (wfc COS 9 ^ v±^ sin ( 



2i;2 



"th.e 



2v' 



"th,e 



(27) 



2wth,i 



4. = + 
cos^? 5 



le 

TOe 



1 1 W 1 

■P— + titP- 



k Vk ' k'^ ' ^ '^'^ 



ITT , 
~k 



2v^ 



2vk 



(vk COS 9 — v± sin ( 



ITTLU dS(vk) 
k^ dvk 

/_ 2«2 \ 



cfc 



2w I eSr 
— H — 



sin 6* 5 



2t;^ 



2w 



(ufe COS 9 — v± sin 6*) 7 



2«^ 



fn.,^ (28) 



for mode (1). Similarly, from equations (|2^) and (P6D, we obtain 

1 „ 1 „ 1 ITT 



fk.t 



le 

TOe 



^ - - ^ - -■■(■/ N iiruj dS(vk) 



V c 



2vfc 



2vth., 







cos 6' 









2vk 



(vk cos — z;_L sin 6') 7 



2i;2 



(wfe cos 9 — v± sin ^ 



2vk 



2wth/ 



2^2 eSr 

2wth,e 

1 „ 1 W „ 1 ITT 
-P \ P 

k Vk k^ v'j. 
2^2 \ 2vk 



th,e 

ie 



2v' 



/m,e; 



(29) 



fc k"^ dvk 



{vk cos ( 



v\ sin (?) 7 - 



2v^ 



"thA 
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ck 



2v, 



[vk COS — sm t^j 7 5 — 



(30) 



for mode (2). 

The Poisson equation is written as 

zfci?fc = 47re((4,)-(4^)). (31) 

Keeping only the first non- vanishing order in e for mode (1), each term on the right-hand side provides 

,2 



2 

-j-Ek- eSr — sin 6* 



2 



pi 



i?z sinf 



(32) 
(33) 



where and Wpi are the electron and the ion plasma frequencies, defined by lu^c = \/ (47moe^/?ne) and 
Ljpi = -\/ (47rnoe^/m,i), respectively; Ad = ■\/A:Br/(47rnoe^) is the Debye length. It is trivial why those 
equations do not contain E± since the transverse component is not constrained from the Poisson equation. 
Note that the first term in the right hand side for ion is exactly the same as that for electron. Therefore, 
the contribution of ion to the charge density is non-negligible. On the other hand, the second term for ion 
is smaller by factor of nie/mi than that for electron. This is the essential point to make possible in-phase 
acoustic oscillation between electron and ion. Then the spatial charge carried by ion acoustic oscillation is 
able to be canceled by electron to keep charge neutrality. In the case of pure electron plasma examined by 
Ramani & Laval (1978), electrons have to keep charge neutrality by themselves. Therefore, the amplitude 
of acoustic oscillation for electron plasma has to be almost zero. The Poisson equation provides the relation 
between electric and magnetic fields as 



eOT 



smf 



Vth.e 



4 l + (fcAD)2 c 



(34) 



The term (fcAn)^ in denominator comes from the left-hand side of the Poisson equation and always very 
small number since the wavelength of the mode which we are interested in is much larger than the Debye 
length. We, therefore, neglect this term. This corresponds to the so-called Plasma approximation (Chen 
1974; Tanaka & Nishikawa 1996) and ensures that the charge neutrality is kept in high accuracy even 
though the longitudinal oscillation exists. When the wave vector is parallel to the temperature gradient 
{9 — Q), this equation tells Ek = and the mode (1) becomes pure transverse. Compared with the case of 
pure electron oscillation as in Ramani & Laval (1978), the right-hand side of equation ( |3^ ) is factor of two 
smaller, which makes mode (1) more unstable as shown in below. No constraint comes from the Poisson 
equation for mode (2) since it is a pure transverse mode. 



We are now in a position to derive the dispersion relation. The Ampere's law is written as 



, „ 47re 
c 



dEj_ 



for mode (1), and 



ikB^ 



Aire 
c 



dE, 
cdt 



(35) 



(36) 
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for mode(2). Under the low-frequeney condition, the last terms in the right-hand side of equations ( pq ) 
and (|3^), those are the displacement current, are negligibly small. We thus neglect these terms. Further, 
the contribution of ion to the current density is order of ^m^jmi smaller than that of electron since the 
current density carried by each particle is proportional to the thermal velocity of each particle. Therefore, 
the current density carried by ion introduces only a negligible contribution to the dispersion relation. The 
dispersion relation for the real parts Wr and the imaginary parts in the leading order of e is then obtained 
as follows: 

<^r = — r"^^th,e cos 0, (37) 



4 



2 



=fct;th,e(2 cos" B - sin^ B)^^[— \ k^vt^^^ (38) 



for mode (1), and 



4^ ' \/7r 



OJr = ^Wth,efcC0s6', (39) 



g2^2 1 / C 1 
= —-^Wth.efc 008^6' -j=[ VtXi.ek^, (40) 



4 

for mode (2). 

Since the imaginary part LOi is of order , we have checked whether the leading order of the dispersion 

("2") ("2") 

relation is changed when the second-order distribution functions /e and fl are taken into consideration. 
We have confirmed that these second-order distribution functions only introduce one order higher terms 
both in real and imaginary parts, and thus main results shown above are not changed. 



5. Instability and mode characteristics 

The dispersion relations (|37|), (|3^), ( p9[ ) and ( ^ ) show that the low- frequency mode, for which the 
phase velocity of the wave is slower than ion thermal velocity, can exist as long as e < ^/m^frni ^ 0.025 is 
satisfied. This is the first confirmation that the Ramani & Laval (1978) type instability can exist beyond 
the wall of square root of the mass ratio. The dispersion relations are almost the same as those obtained for 
electron plasma by Ramani & Laval (1978). However, there is a slight difference in imaginary part of the 
mode (1). The difference comes from the difference in the charge neutrality condition (Q). As explained 
in §4, in our case, the in-phase acoustic oscillation between electron and ion is possible. Hence, the spatial 
charge carried by ion acoustic oscillation is able to be canceled by electron to keep charge neutrality. 
However, in the case of pure electron plasma examined by Ramani & Laval (1978), electrons have to keep 
charge neutrality by themselves. Therefore, the amplitude of acoustic oscillation for electron plasma has to 
be almost zero. 

The characteristics of the instability are summarized as follows. For mode (1), the imaginary part of 
the wave frequency is positive and the instability sets in, when the direction of the wave vector is within the 
double cone spanned by G [— ^'cr, ^cr] and e [tt — 0cr, ti" + ^cr], where 0^ = arccos(l/-\/3) (0 < 6'cr < ''■/2). 
For comparison, we re-calculated the dispersion relations for pure electron plasma. In the case of pure 
electron plasma, a factor of 2 appears in front of the sin" B term in the imaginary part (this result is 
slightly different from Ramani & Laval (1978)'s result.). Therefore, the unstable region in /c-space in the 
case of pure electron plasma is somewhat narrower than our case. It shows that acoustic oscillation of ion 
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is assisting the instability. For mode (2), the unstable mode exists for all direction of the; wave vector, 
although the growth rate decreases as increase from to 7r/2 for fixed k. Two modes get identical when 
6 — 0. For both modes, the growth rate is maximum when = 0. Since the imaginary part of the wave 
frequency is the third order polynomial of k with negative coefBcient for k^, there exists a maximum growth 
rate Wi^max- For mode (1), Wi^max is given as 



e^'^T Vth,e _ |Q„_2fl i|3/2 

12\/37f c 



Wi,max ~ ^ —f-LOpe |3cOS - 1| , (41) 



when 

k = k^,^ = ^'^\3cos'0-l\'^\ (42) 
2V3 c ' ' 



for 6 e [-^cr, ^cr] or ^ e [tt - ^cr, TT + ^cr] • For modc (2), 

i,e 



e^6^ Vth e I 3 /)| rAO\ 
Wi,max ^ T-7= -l^pe CDS'" 0\ (43) 



when ^ 

k = fcmax = — ^ — |C0S 6I| , (44) 

for arbitrary 6. For the wave vector which satisfies the above conditions, the real part cji of the wave 
frequency is obtained as 

~ ^ —IjJpe 1 3 cos y — 1 1 cos (45) 

8v 3 c 

for mode (1), and 

1 2 r2 ^th,e 2 /) 

Wi ;=£ ^WoeCOS 6' (46) 

4\/6 c 

for mode (2). 

The electric field strength can be related to the magnetic field strength. For mode (1), 

l^^fcl^^ ^^l-Bfcllsm^l. (47) 

When k = /Cmax, the electric field perpendicular to the /e-direction is related to the magnetic field as 

I^^|-^^l4ll-s^|. (48) 

For mode (2), 

\E^\^'-^!^\B,\\cos0\, (49) 



The electric field strength is then order of e smaller than the magnetic field strength. This nature 
combined with the low-frequency nature of the mode shows that the mode has similarity with the 
magneto-hydrodynamical mode. However, the magneto-hydrodynamical treatment cannot identify the 
mode and the instability. For example, the instability is microscopical, in which case the resonance of 
particles with waves is essential. 
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The possible reduction factor of the mean free path, in other words the possible lowest value of e, is 
set by the application limits and assumptions summarized in §2. In the following discussion, we treat the 
fastest growing mode, that is 9 = case. First, the condition (^) provides 

1 C \ ' / 



Lcupcj V 1.0 



^>4(t— ) (t^) • (50) 

Second, the age condition provides 



C7 \ 1/3 / 



e>2 — — {-^ . (51) 



Third, the coUisionless condition is satisfied when 27r/fc,jiax < e x L, which yields 

-, \ 1/2 / c X -1/2 

1 c \ I Ot ^ 



Finally, the condition that the life time of the interface given by equation (|l|) must be longer than the age 
of the cool region, iagc(cool), provides the upper limit on e as 

-1 / T \ / 1 rr^ \ -0.5 

(53) 



' IQiOyr J \AkpcJ \AkeY J 



6. Application to the cluster of galaxies 

6.1. A jam bun model 

For the several nearby clusters of galaxies, the existence of at least the two different temperature 
gas in the central region is spectro-scopically confirmed (Ikebe 1995). If the cool gas found in the cluster 
central region is a single block of the low temperature gas located at the bottom of the cluster gravitational 
potential well like a jam in a jam bun, the interface between the cool and hot gas can be well described by 
the pressure equilibrium condition and the results obtained in this paper can be applied. One of the well 
studied such clusters is the Centaurus cluster (Ikebe et al. 1999). In the following discussion, the cool gas 
found in the Centaurus cluster is assumed to be well described by the jam bun model. The temperature of 
the cool region is about fceT'i ^ IkeV and that of the hot region is about fcBT'2 ^4keV in the Centaurus 
cluster (Ikebe et al. 1999). The factor St is obtained as 5t ~ (T2 - Ti)/(T2 + Ti)/2 ~ 1.2. The size 
of the region, where the cool component exists, is about the size of the X-ray halo around the elliptical 
galaxies, say 50-100/i^^kpc (Forman, Junes, & Tucker 1985; Matsushita et al. 1998; Ikebe et al. 1999). 
Using the best fit values for the electron density profile for the hot phase gas obtained by ASCA results 
(Ikebe et al. 1999), the electron number density of the hot phase gas at around lOOkpc can be estimated as 
2.5 X lO^^cm^'^. Then the electron mean free path Ac due to Coulomb collisions is about 2kpc. Hence, the 
condition (^ implies L > 2kpc, at least when the low-temperature region was formed. The age condition 
(|l|) leads to 



The coUisionless condition (p3) also provides a lower limit on e as 
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Therefore, the age condition yields the lower limit on e when L > 0.25kpc and is compatible with the life 
time condition (|53|). On the other hand, the condition jssl ) provides the lower limit on e when L < 0.25kpc. 
In this case, for this condition to be compatible with the life time condition (|53|), the lower limit of the size 
of the interface region is set as 

I 2.5 X 10-3cm-3 j (ik^j (12; "^P" 

Compared with the cluster pressure scale height, the obtained small values of the lower limit on L both from 
the age condition and the condition (|5^) guarantee the application of the pressure equilibrium condition. 
Therefore, we conclude that the instability found in this paper can be a possible mechanism to inhibit the 
heat conduction between cool and hot phases in the cluster central region, if the cool gas is a single block 
of the low temperature gas located at the bottom of the cluster gravitational potential well. 



6.2. The raisin bread model 

If the cool gas is contained in small blobs and the blobs distribute in the entire cluster central region 
like raisins in a raisin bread, the obtained results cannot be simply applied. The cooling flow model, which 
is one of the widely- accepted model to describe the cluster central region (Fabian 1994), assumes that the 
ICM in the cluster central region is in multi-phase and is composed of many small cool clumps (raisins) 
which have different temperatures from clump to clump. Even in the two temperature phase model (Ikebe 
et al. 1999), the raisin bread model is thought to be more plausible to describe the situation in the cluster 
central region. It is well known that these cool clumps, which is in pressure equilibrium with surrounding 
hot medium in the cluster central region, oscillates owing to buoyancy force (Balbus 1986; Malagoli, Rosner 
and Bodo 1987; Balbus and Soker 1989) and that they are fragmented (Nulsen 1986; Balbus 1986) and 
merge into the surrounding hot medium within a few buoyancy oscillation period, that is ~ a few x lO^yr 
(Hattori and Habe 1990; Malagoli, Rosner and Fryxell 1990; Reale et al. 1991; Yoshida, Hattori and Habe 
1991). It has been shown that the magnetic field frozen in the ICM can prevent the buoyancy oscillation of 
the cool clumps (Nulsen 1986; Loewenstein 1989; Balbus 1991; Hattori, Yoshida and Habe 1995). Therefore, 
the raisin bread model requires the existence of the magnetic field in the plasma. Nulsen (1986) pointed out 
that if the Alfven velocity va of the hot medium is larger than the terminal velocity vt of the cool clumps, 
that is, if 

VA>vt, (57) 

the relative motion of cool clumps against the surrounding medium is stopped and the fragmentation of the 
clumps are prevented. The terminal velocity of a cool clump with a relative density excess of 5p and a size 
of r is given by vt = ^ {Sp/p){r/R), where R is the pressure scale height of the cluster. It was confirmed 
by the numerical simulation (Hattori, Yoshida and Habe 1995) that this condition works well even for the 
cool clumps with 6p/ p > 1. It is important to note that even if the magnetic field line is parallel to the 
gravitational force, the buoyancy oscillation can be prevented by the magnetic tention which is first pointed 
out by Loewenstein (1989). 

It can be shown that the gas pressure equilibrium condition is a good approximation even for the 
cool clumps in the raisin bread model. Suppose that a cool gas clump is in pressure equilibrium with 
surrounding medium which is immersed in the background magnetic field. The cool gas clump falls ahead 
of the surrounding medium owing to its excess weight. As a result, the magnetic field configuration of the 
surrounding medium is locally disturbed. The infall of the cool clump is stopped, when the force produced 
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by the disturbed magnetic field acting on the cool clump balance with the buoyancy force. This condition 
is written as 

(58) 



r ^ " R 



where Pm and Pg are magnetic and gas pressures of the hot surrounding medium, respectively; 5b = SB/B 
and dp — Sp/ p are deviations of the magnetic field strength and the gas density in the clump, respectively, 
relative to those of surrounding medium when the force balance is realized. This condition can be rewritten 
using the plasma /?(/? = Pg/P,n) as 

SB^f3Sp^. (59) 



From equation (|59|), the condition (57) can be interpreted as 1 > Sb, and the condition for the possible cool 
clumps in the fixed-/? plasma is obtained as 

In the equilibrium state, the total pressure in the clump must be the same as the total pressure of the 
surrounding hot medium. This condition provides 

\SPg\ ~ \6BPml (61) 

in the clump. Therefore, the deviation of the gas pressure relative to the surrounding gas pressure can be 
estimated as 

\SPJ , r 



ai 



Pg '^R- (^2) 

From equations ( |60| ) and (^), we obtain the constraint on the relative deviation of the gas pressure: 

Since the hot gas in the cluster is high-/? plasma, say /? ^ 10 — 100 (Kim et al. 1989; Makino 1997), the 
relative deviation of the gas pressure is at most 10% and very small. 

The arguments presented above ensures that the model of the interface adopted in this paper is 
applicable to the case when the raisin bread model is the good representative for the cluster central region. 
Therefore, the results obtained in this paper can provide some hint for reduction of the heat conductivity in 
cooling flow. However, the further analysis including the background magnetic field is required to apply the 
raisin bread model since the background magnetic field introduces another resonance related to cyclotron 
frequencies of electron and ion. 



7. Discussion 

We have shown that there exists low-frequency growing modes driven by a global temperature gradient 
in electron and ion plasmas. The instability provides a hint for reduction mechanism of the thermal 
conductivity to explain the observational evidence of co-existence of the multi-phase gas in clusters of 
galaxies. However, to answer whether the instability really reduces the heat conductivity down to the 
required level, non-linear saturation level of the instability and the effect of the back ground magnetic field 
must be examined. Although Levinson & Eichler (1992) has studied the effect of the magnetic field on the 
Ramani and Laval type instability for pure electron plasma, there is no study which includes ion with back 
ground magnetic field. 
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The evolutionary nature of the modes in the interface region can be briefly described as follows. When 
the low-temperature gas is injected into the intracluster medium, any wave is not excited in the interface 
region yet. Hence, e is defined by the mean free path due to Coulomb collisions and is order of 1 to 0.1, as 
shown in §1. At this initial time period, the Ramani & Laval (1978) type instability excites the waves only 
for electron plasma, since the phase velocity of excited waves is about the electron thermal velocity and thus 
much faster than the ion thermal velocity. As waves grow, the mean free paths of electron and ion become 
shorter and e gets smaller. When e goes down to ^Jm^Jnii^ ions start to respond the waves, and the modes 
switch to those described in this paper. Since the growth rate of the instability is also proportional to 
and very large at the initial stage, the condition ( |50|) is satisfied. We thus do not have to worry about the 
disappearance of the interface due to the thermal conduction before the waves grow. The waves must be 
scattered back at the edge of the interface and must be confined in the interface region, because the waves 
driven by the instability cannot propagate in the uniform temperature region. 

Since the instability creates and grows the magnetic field without any seed field, it provides a new 
possibility for the origin of the cosmic magnetic field. The Biermann battery effect (Biermann 1950) is the 
mechanism which has been considered as one of the central mechanism to create the seed magnetic field in 
our universe (Kulsrud et al. 1997). The Biermann battery effect requires that the direction of the density 
gradient of the gas is not parallel to the direction of the temperature gradient. However, the instability 
found in this paper requires only the existence of the temperature gradient. In the astrophysical situation, 
the situation that two different temperature plasmas are contacting keeping pressure equilibrium can be 
commonly found compared with the situation that the Biermann battery effect may act. Therefore it may 
be interesting to examine how the instability found in this paper plays a role on the origin of the cosmic 
magnetic field. 

Finally, it is interesting to note that our results suggest that the realization of the global thermal 
equilibrium is postponed by the local instability which is induced for quicker realization of local thermal 
equilibrium state in plasmas. This result is of course not trivial. Once the global thermal equilibrium 
is achieved, the local thermal equilibrium is realized simultaneously. Therefore, letting electrons be free 
to conduct thermal energy without any instability is the fastest way to achieve thermal equilibrium. 
However, the plasma prefers restoring the local deviation of the velocity distribution function from the 
Maxwell-Boltzmann and excites waves by the instability. As a result, the realization of the local thermal 
equilibrium is also postponed. The plasma be hurried into an error. 

Authors would like to thank Masahiro Hoshino, Rudolf Treumann, Yousuke Itoh, Yutaka Fujita, 
Masahide lijima for many instructive and fruitful discussion on this topics. 

REFERENCES 

Balbus,S.A. 1986, ApJ, 303, L79 
Balbus,S.A., & Soker,N. 1989, ApJ, 341, 611 
Balbus,S.A. 1991, ApJ, 372, 25 
Biermann, V. L. 1950, Z. Naturforsch, 5a, 65 

Chapman, S., & Cowling, T. G. 1960, The Mathematical theory of non- uniform gases, (Cambridge 
University Press, Cambridge) 

Chen, F. F. 1974, Introduction to Plasma Physics, (Plenum Press, New York) 



- 15 - 



Fabian, A. C. 1994, Ann.Rev.Astron.Astrphys., 32, 277 
Forman, W., Jones, C, & Tucker, W. 1985, ApJ, 293, 102 
Fried, B. D. 1959, Phys. Fuild., 2, 337 

Fukazawa, Y., Ohashi, T., Fabian, A. C, Canizares, C.R., Ikebe, Y., Makishima, K., Mushotzky, R. F., & 

Yamashita, K. 1994 PASJL 46, L55 

Hattori,M., & Habe,A. 1990, MNRAS, 242, 399 
Hattori, M., Yoshida, T., & Habc, A. 1995, MNRAS, 275, 1195 
Hattori, M., Kncib, J.-P., & Makino, N. 1999, PTP Suppl., 133, 1 
Ikebe, Y. 1995, PhD. Thesis, (University of Tokyo) 

Ikebe, Y., Makishima, K., Fukazawa, Y., Taniura, T., Xu, H., Ohashi, T., & Matsushita, K. 1999, preprint 

Kim,K.T., Kronbcrg,P.P., Giovannini,G., & Vcnturi,T. 1989, Nat, 341, 720 

Kulsrud, R. M., Ccn, R., Ostriker, J. P., & Ryu, D. 1997, ApJ, 480, 481 

Levinson, A., & Eichler, D. 1992, ApJ, 387, 212 

Locwcnstcin, M. 1990, ApJ, 349, 471 

Makino,N. 1997, ApJ, 490, 642 

Makishima, K. 1997, in Imaging and Spectroscopy of Cosmic Hot Plasmas, eds. F. Makino & K.Mitsuda 
(University Academy Press, Tokyo), pl37 

MalagoH,A., Rosner,R., & Bodo,G. 1987, ApJ, 319, 632 

MalagoH,A., Rosner,R., & Fryxell,B. 1990, MNRAS, 247, 367 

Matsushita, K., Makishima, K., Ikebe, Y., Rokutanda, E., Yamasaki, Y. N., & Ohashi, T. 1998, ApJ, 499, 
L13 

Melrose, D. B. 1986, Instabilities in space and laboratory plasmas, (Cambridge University Press, Cambridge) 

Nulscn, P. J. E. 1986, MNRAS, 221, 377 

Pistnner, S., Levinson, A., & Eichler, D. 1996, ApJ, 467, 162 

Ramani, A., & Laval, G. 1978, Phys. Fluid., 21(6), 980 

Reale,F., Rosncr,R., MalagoU,A., Pcrcs,G., Serio,S. 1991, MNRAS, 251, 379 

Sarazin, C.L. 1988, X-ray emission from clusters of galaxies, (Cambridge University Press, Cambridge) 
Sokcr, N., & Sarazin, C.L. 1990, ApJ, 348, 73 

Tanaka, Y., Inouc, H., & Hoh, S. S. 1994, Publ.Astron.Soc.Jpn., 46, L37 

Tanaka. M., & Nishikawa, K. 1996, Physics of high temperature plasmas, (Maruzen, Tokyo) 

Wcibcl, E. S. 1959, Phys.Rev.Lctt., 2, 83 

White, R. A., & Sarazin, C. L. 1987, ApJ, 318, 612 

Yoshida,T., Hattori,M., Habe,A. 1991, MNRAS, 248, 630 



This preprint was prepared with the AAS lAT^jX macros v4.0. 



- 16 - 



Fig. 1. — A sketch of one-dimensional distribution of temperature T[x) in the interface region. Ti and T2 are 
the temperatures in the cool and hot regions, respectively. Also shown are the number density distribution 
n{x) (dashed line) and the distribution of pressure (dotted line), which is constant along the a;-direction. 



